Melting of alloys along the inter-phase boundaries in eutectic and peritectic systems 
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We discuss a simple model of the melting kinetics along the solid-solid interface in eutectic and 
peritectic systems. The process is controlled by the diffusion inside the liquid phase and the existence 
of a triple junction is crucial for the velocity selection problem. Using the lubrication approximation 
for the diffusion field in the liquid phase we obtain scaling results for the steady-state velocity of 
the moving pattern depending on the overheating above the equilibrium temperature and on the 
material parameters of the system, including the dependences on the angles at the triple junction. 



I. INTRODUCTION 

The systematic investigation of melting kinetics in al- 
loys, and particular in eutectic and peritectic systems, is 
much less developed than the investigation of solidifica- 
tion (for a recent review on solidification see, for example 
PJ and references therein). Microstructures, being at the 
center of materials science and engineering, are formed 
during the solidification process and, in this sense, the 
melting process is less attractive for practical applica- 
tions. 

However, the interfacial pattern selection problem dur- 
ing the melting process might be very interesting. For 
example, in our previous publications 0, Q we discussed 
a free boundary problem for two moving interfaces that 
strongly interact via the diffusion field in the liquid layer 
between them. This problem arises in the context of liq- 
uid film migration during the partial melting of solid al- 
loys and could also be relevant to the sintering process 
in the presence of the liquid phase For the melting 
of one-phase alloys to proceed in this way, the local equi- 
librium concentrations have to be different for the two 
interfaces providing the driving force for theprocess. ft 
is by now well accepted (see, for example, @,Ll) that the 
difference of the equilibrium states at the melting and 
solidification fronts is due to the coherency strain energy 
which is important only at the melting front because of 
the sharp concentration profile ahead the moving melting 
front. 

The other source of elastic deformations during the 
melting process, even in pure materials, is the density 
difference between the solid and the liquid phase. If the 
melt inclusion is entirely inside the solid matrix, inhomo- 
geneous elastic deformations inevitably arise. The pecu- 
liar behavior of the melting kinetics in such systems was 
discussed in || . 

The main purpose of this paper is to describe the prob- 
lem of contact melting in eutectic and peritectic systems 
along the boundary between two solid phases (see Fig. 1 
and Fig. 2). The local concentrations at the L/a and 
L/P interfaces in such systems are different because of 
the chemical difference between the a and (3 phases, and 
weak coherency strain effects are not so important here. 
If we also assume that the liquid phase extends up to 
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FIG. 1: Schematic presentation of the phase diagram a) and 
configuration of different interfaces near the triple junction 
b) in eutectic systems. x a p(y) is the interface between two 
solid phases ct and (3; the interfaces between the liquid phase 
L and the solid phases are denoted by x a (y) and Xp(y). In 
the steady state regime this configuration moves along the y 
axis with a constant velocity v. The origin of the coordinate 
system is located at the triple junction. 



the sample surfaces, the mentioned elastic deformations 
due the density difference are not important either, be- 
cause a weak hydrodynamic flow inside the liquid phase 
compensates the density difference. 

We concentrate here on the velocity selection prob- 
lem during the melting along the solid-solid interface. 
The presence of the triple junction (see Fig. lb and Fig. 
2b) plays a crucial role in this process. In the classical 
problem of dendritic growth, where the triple junction is 
not present, the velocity selection is controlled by tiny 
singular effects of the anisotropy of the surface energy 
(for review see 0, Hlj). The triple junction produces 
a very strong perturbation of the liquid-solid interfaces 
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and controls the velocity selection. The other important 
difference compared to the classical dendritic problem is 
that the kinetics of the contact melting in eutectic and 
peritectic systems is controlled by the diffusion inside the 
needle-like liquid phase and not in the outer phase. 

All features of the contact melting process mentioned 
so far are common for both systems. There is, however, 
important difference between eutectic and peritectic sys- 
tems. In the eutectic system (Fig. 1) both interfaces 
L/a and L/(3 are melting fronts while in the peritectic 
system (Fig. 2) the L/a interface is a solidification front 
if the temperature To is above the peritectic tempera- 
ture Tp. In other words, during the melting of peritec- 
tic systems the low-temperature /3-phase melts while the 
high-temperature a-phase solidifies. It means that, ad- 
ditionally to the formation of the liquid phase, the poly- 
morphic transition (3 — -> a occurs. In this context the 
contact melting in the peritectic system is similar to the 
process of liquid film migration, mentioned above, where 
the liquid film is also located between melting and solid- 
ification fronts. 

Finally, we discuss in this paper the evolution of the 
solid-solid interface due to surface diffusion. This process 
inevitably arises due to the deviation of the interface from 
a fiat configuration in the vicinity of a triple junction. 
We show, however, that surface diffusion does not play a 
controlling role in the melting kinetics, and only allows 
for the necessarily adjustment of the solid-solid interface. 



II. ISOTHERMAL MELTING IN EUTECTIC 
SYSTEM 

We consider the two-dimensional problem of the simul- 
taneous melting of two eutectic phases a and (3 along the 
boundary between them. The phase diagram and the 
configuration of different interfaces near the triple junc- 
tion are schematically presented in Fig.l. The compo- 
sitions of the solid phases are close to their equilibrium 
values C a and Cp at the eutectic temperature. We as- 
sume that diffusion in the solid phases is very slow. The 
temperature of the sample, To, is slightly above the eu- 
tectic temperature Te (for the notations see Fig.l). At 
small overheating the concentration field in the liquid, 
C(x,y) obeys the Laplace equation 

AC = 0. (1) 

We assume local equilibrium at the liquid-solid interfaces. 
Then the concentrations at the L/a and L//3 interfaces 
are: 

C La {y)=Ci a + {Cp-C a )d a x'^y), 
C L p(y) = C e L p + (Cp-C a )dpx'p , (y). (2) 

Here the notations of the different concentrations are 
clear from the phase diagram (Fig. la); the capillary 
lengths are 
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where c a (,8) is the surface energy of the L/a(f3) interface, 
O is atomic volume which is assumed to be the same in all 
three phases, /l(C) is the free energy of the liquid phase 
per atom and f'l(C) is its second derivative with respect 
to C. The mass balance conditions at these interfaces 
read: 



DdC/dx = v(CE-C a )x' a (y), 
DdC/dx = -v(Cp-C E )x' (y), 



(3) 



where D is the diffusion coefficient in the liquid phase. 

Eqs. ©-© ar e written for the case of small angles 
9, ip a and tpp at the triple junction. This leads to the 
small angles along the whole interface, x a -C 1 and x^ <C 
1. In this case we can use the so-called "lubrication" 
approximation for the solution of the Laplace equation. 
In this approximation one neglects the derivatives with 
respect to the "slow" y variable as compared with the 
derivatives with respect to the "fast" x variable: 



C(x,y) = C (y)+B(y)x. 



(4) 



The slow variable functions Co(y) and B(y) can be 
found from the boundary condition, Eq. J3J): 



Co (y) 



C La (y)xp(y) - C L p(y)x a (y) 
xp(y) - x a (y) 

B (y) = c ^(y)-cUy)^ (5) 



xp{y) - x a (y) 
From Eqs. 10 and (0J we find the following relations 



, _ D B(y) 

Ja v{c E -c a y 



, D B(y) 

13 ~ v(c p -c E y [ ) 



which together with Eqs. (JHJ and form a closed sys- 
tem of two second order differential equations for the 
front profiles x a (y) and Xp(y). They are subject to the 
boundary conditions at the triple junction: 



x a {O) = x (O) = O, 



(7) 



x a (0) = ip a - 6 = a ; x f3 (O) = -(( Pl3 + 6) = -6 (8) 

(see Fig. lb). It follows from Eqs. JBJ and Q that 

x a (y) = -gxp(y), (9) 

where g = (Cp — Ce)/(Ce — C a ). We also find from Eq. 
© that 



1 + 5 



1+3 



1 + 3 



\Ce - C a (p)\(Cp - C a )f'l(C E )' 



(10) 

The relation 10 allows to eliminate the profile x a and to 
write the closed equation for the profile xp. Moreover, 
one can integrate this equation once with respect to y and 
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finally obtains a nonlinear first order differential equation 
for the profile xp: 



(dp + gd a )x /3 = (dp 



gd a )9p 
C a 



Ay 



vg 



i . 



(11) 



2D 

is the dimcnsionless 



Here A = (C% p - Cl a )/(Cp_ 
overheating above the eutectic temperature. 
Let us introduce the dimcnsionless coordinates 



X 







xpA 



Y = 



yA 



'p\u>p -r gd a ) 
Then Eq. (fTT|) takes the form 



Op (dp + gd a ) 



X 







1 + Y 



uX} 



with the dimcnsionless parameter v: 

g(dp + gd a )9l 

v = v m • 

2DA 2 



(12) 



(13) 



(14) 



The solution of this equation starts at the origin with 
Xp(0) = —1 and should have parabolic asymptotics 

Xp = \J\Y\jv for large values of \Y\. It turns out 
that such a smooth solution exists only if the parame- 
ter v = v* w 1.06. In other words, this is a nonlinear 
eigenvalue problem which leads to the selection of the 
velocity v: 



2DA 2 



g(dp + gd a )9l 



(15) 



We see that the posed problem has a relatively simple, 
essentially analytical solution. 



III. ISOTHERMAL MELTING IN PERITECTIC 
SYSTEM 

The phase diagram and the configuration of different 
interfaces near the triple junction are schematically pre- 
sented in Fig. 2. The temperature of the sample is To and 
slightly above the peritectic temperature Tp (for the no- 
tations see Fig. 2a). Overheating is assumed to be small 
and the compositions of the solid phases are close to the 
equilibrium values C a and Cp at the peritectic tempera- 
ture. Eqs. C0l-l|IU are valid also for the peritectic system 
where Ce should be replace by Cp and in the second 
equation in @ one should change the sign in front of the 
capillary term. 

C L p(y) = Cip-(Cp~C a )dpx"p(y). 

This reflects the fact that the equilibrium between the 
/3-phasc and the liquid phase corresponds to the negative 
slope of the liquidus line for the peritectic phase diagram 
and to the positive slope for the eutectic diagram. As a 
result, Eqs. l|9")l- lllU|) remain the same and Eq. Hlfl is 
replaced by 



(dp - gd a )xp = (dp - gd a )9p - Ay 



yg 

2D 



0' 



(16) 



a) 




b) 





y 


a 1 


P 


\ 


V 




X 








— \'e x 




X M l N 



FIG. 2: Schematic presentation of the phase diagram a) and 
configuration of different interfaces near the triple junction 
b) in peritectic systems. x a p(y) is the interface between two 
solid phases a and /3; the interfaces between the liquid phase 
L and the solid phases are denoted by x a (y) and xp(y). In 
the steady state regime this configuration moves along the y 
axis with a constant velocity v. The origin of the coordinate 
system is located at the triple junction. 



We note that the parameters g and A for the peritectic 
system (Fig. 2) are negative if we define them in the same 
way as for the eutectic system: 

g = (Cp - Cp)/(Cp - C a ), A = (Clp - CU/(Cp - C a ) 

Thus, replacing — g by \g\ and —A by |A|, we recover Eq. 
(|T3j) with 



X, 



xp\A\ 



\(dp + \g\d a 



Y 



W|A| 



9p(dp + \g\d a ) 



and 



\g\(dp + \g\d a )9 



<p 



1 



I. '/I 



2L>A 2 ' 
Finally, for the steady-state velocity we obtain 

2D A 2 



\g\(dp + \g\d a )9l 



(17) 



(18) 



(19) 



As we have already noted in the Introduction, the in- 
terface x a = \g\xp represents now a solidification front. 
During the melting, in the peritectic system, the low- 
temperature /3-phase melts while the high-tcmpcraturc 
a-phasc solidifies. In addition to the formation of the liq- 
uid phase the polymorphic transition /? — > a takes place. 
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IV. EVOLUTION OF THE SOLID-SOLID 
INTERFACE 

We have already mentioned that diffusion along the 
solid-solid interface plays an important role allowing the 
necessary rotation of the structure in the vicinity of the 
triple junction, x' a g(0) = —9 (see Fig. lb and Fig. 2b). 
At the same time, far away from the triple junction, the 
solid-solid interface is parallel to the y axis (x' a p(oo) = 0). 

The corresponding surface diffusion equation which de- 
scribes the x a p(y) interface has been discussed in [l2T | 
and, essentially, it is given by the classical Mullins equa- 
tion [13: 

2d a pd a/ 3D a p—j-f- - v ~^- = I 20 ) 

In this description the diffusion takes place inside of layer 
of the thickness 25 a p (diffusion coefficient is 2? a( a) and 
d a s is the corresponding capillary length: 

a0 2{Cp-c a yfi>{c a )f»{Cpy 

Here f a (/3)(C) is the free energy of the a((3) phase per 
atom. The solution of this equation which has a proper 
behavior at large values of y and satisfies the boundary 
conditions x a p(0) = 0, x aj3 {Q) = —9 and x a/3 (0) = k(0), 
is 

x a p(y) = A ~ cx P(~yJ [^cos(— — ) +Bsin(— )], 

A=[ K (0)/q-6]/q; B = [k(0)/ q + 9\j (q\/~3), (21) 

where q 3 = v / {25 a pd a pD a p). This solution is uniquely 
defined because the parameters v, 9 and k(0) are known 
from the solution given in the previous sections. The 
velocity is given by Eqs. (|15|) and (|19fl for the eutectic 
and peritectic systems respectively; the angle 9 is given 
by Eq. (|10f> for both cases and the curvature of the a/f3 
interface at the triple junction, n(0) ~ A/d a p, is also 
known since we have already found the concentrations at 
the triple junction. The procedure of calculation of k(0) 
is straightforward but tedious and involves new thermo- 
dynamical parameters. Therefore, we explain this proce- 
dure only schematically: The concentrations in the solid 
phases at the triple junction can be calculated, on one 
hand, using the condition of local equilibrium with the 
liquid phase (along the x a ^ interface) and, on the other 
hand, using the condition of local equilibrium between 
two solid phases (along the x a p interface). At the triple 
junction these compositions should coincide. The capil- 
lary corrections to the equilibrium solid concentrations 



can be written similar to Eq. Q for all there interfaces. 
Since the profiles x a and xp have been already found, 
the mentioned continuity condition allows to calculate 
the curvature k(0) of the solid-solid interface. 

We note, however, that the term ft(0)/g in Eq. (|21|l is 
of order 9{D aP A/D) 1 / 3 and can be neglected compared 
to 9 for small A. This means that the tedious calculations 
explained above are not needed in the limit of small A. 
We also note that the characteristic length scales of the 
solid-solid intcrfacial pattern are small compared to the 
characteristic length scales of the melt structures by the 
same small parameter (D a pA/ 'D) 1 / 3 . Moreover, the dif- 
fusional flux along the solid-solid interface has a nonzero 
value at the triple junction. In principal, this flux has to 
be taken into account in the description of the diffusional 
field in the liquid phase. This effect has been neglected in 
the previous sections. More careful analysis shows that 
corrections to Eq. (6) due to this effects are small by 
the same small parameter, (D a pA/ 'D) 1 / 3 « 1. Thus, in 
the limit of small A the surface diffusion process has no 
influence on the kinetics of the contact melting and al- 
lows only for the necessarily adjustment of the solid-solid 
interface inside the relatively small region in the vicinity 
of the triple junction. 



V. CONCLUSION 

We have developed and analyzed a relatively simple 
model for the melting kinetics along the solid-solid inter- 
phase boundaries in eutectic and peritectic systems. The 
process is controlled by the diffusion inside the liquid 
phase and the existence of the triple junction is crucial 
for the velocity selection problem. The additional as- 
sumption of small opening angles tp at the triple junction 
plays only a technical role and allows to solve the posed 
problem essentially analytically using the lubrication ap- 
proximation for the diffusion field in the liquid phase. 
The obtained scaling results for the dependence of the 
steady-state velocity of the moving pattern on the over- 
heating, which are exact in the limit of small angles, can 
still be used for the moderate values of the angles with 
only prefactors of the order of unity missing. We note the 
obtained scaling relation, v ~ DA 2 /d, is similar to the 
scaling relation in eutectic growth. Using Eq. I|15fl an d 
characteristic values of the parameters D ~ 10 _9 m 2 /s, 
d ~ 10~ 9 m, A ~ 10 -3 and 9 ~ 1 we estimate the char- 
acteristic velocity to be of the order of v ~ 10 -6 m/s. 
We hope that our results will stimulate some new model 
experiments on contact melting phenomena. 

We acknowledge the support by the Deutsche 
Forschungsgcmcinschaft under project SPP 1120. 
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